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Abstract 

In-medium sum rules following from the chiral charge algebra of QCD are reviewed, 
and new sum rules are derived. The new sum rules relate the I G {J PC ) = 1 _ (0 ++ ) 
excitations (quantum numbers of ao(980)) to the scalar and isovector densities, 
and are nontrivial for the isospin-asymmetric medium. We present an extensive 
illustration of the sum rules with help of quark matter in the Nambu-Jona-Lasinio 
model. Collective excitations different from the usual meson branches (spin-isospin 
sound modes) are shown to contribute significantly to the sum rules and to play a 
crucial role in the limit of vanishing current quark masses. 
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1 Introduction 



Over the past years intense efforts have been made to better understand the proper- 
ties of nuclear systems under extreme conditions [1,2]. It is commonly accepted that 
basic properties of hadrons undergo severe modifications in nuclear medium [3-10] We 
expect that at sufficiently large densities chiral symmetry is restored. Moreover, we 
know that already at nuclear saturation density we should find strong medium effects. 
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For instance the quark condensate (qq) is estimated to drop to about 70% of its vac- 
uum value at the nuclear saturation density, as follows from the model-independent 
prediction of Refs. [11,12]. The change in this basic scale of strong interactions, as 
well as other matter-induced effects, undoubtedly lead to severe modifications of in- 
medium hadron properties, whose excitation energies, widths, coupling constants, 
size parameters, etc. undergo changes. The experimental evidence for these effects 
can be found in studies of mesonic atoms, or in the measurements of dilepton spectra 
in heavy-ion collisions in the Ceres [13] and Helios [14] experiments at CERN. Much 
more accurate data on hot and dense matter will be provided by the Hades exper- 
iment, and by RHIC in the near future. It is therefore an important task to better 
understand and describe theoretically mesonic excitations in dense and hot systems. 

Recent years have brought new interesting ideas and developments in this field. The 
incomplete list, relevant for the subject of this paper, contains the possibility of S-wave 
kaon condensation in nuclear matter [15-18], and the application of chiral effective 
Lagrangians and models [18-26] to nuclear systems. General mo del- independent pre- 
dictions for excitations with quantum numbers of the pion, based on chiral charge 
algebra, were made in Refs. [18,27-30]. Our present work summarizes and further 
extends the results presented there. 

The purpose of this paper is twofold: In the first part we review the previously- derived 
current-algebraic sum rules for pionic excitations in nuclear medium (the generaliza- 
tion of the Gell-Mann-Oakes-Renner relation [18,27], the sum rule of Ref. [28]), as 
well as derive new sum rules concerning the excitations with quantum numbers of 
the 5 (or a ) meson (I G (J PC ) = 1~(0 ++ )) (Sec. 2). We discuss formal predictions 
following from these sum rules (Sec. 3). Particular attention is drawn to nuclear mat- 
ter with isospin asymmetry, since this is the case where nontrivial conclusions can 
be drawn for the behavior of mesonic excitations in the limit of vanishing current 
quark masses. We discuss the appearance of very soft modes in this limit. In the pion 
channel there exists a positive-charge mode (for medium of negative isospin density) 
whose excitation energy scales in the chiral limit as the current quark mass itself, 
and the the square root of it, as is the case of the vacuum. In the S channel there 
exists a positive-charge mode (for medium of negative isospin density) whose exci- 
tation energy scales as the difference of the current masses of the u and d quarks. 
These modes are shown to completely saturate the sum rules in the limit of vanishing 
current quark masses. 

In the second part of the paper (Sec. 4-8) we present an extensive illustration of the 
general results with help of quark matter in the Nambu-Jona-Lasinio model [31]. 
Although quark matter is not a realistic approximation to nuclear matter (except, 
perhaps, at very large densities), the model is good for the present purpose. The 
reason is that the Nambu-Jona-Lasinio model is consistent with chiral symmetry and 
complies to chiral charge algebra relations leading to the sum rules. We show that 
the results of the model are highly nontrivial: collective states appear in isospin- 
asymmetric medium (spin-isospin sound modes) and these states are necessary to 
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saturate the sum rules. For certain choice of model parameters, these sound modes 
become the very soft modes in the limit of vanishing current quark masses, and they 
completely saturate the sum rules. Finally, we remark that the Nambu-Jona-Lasinio 
model is interesting in its own, and that much of the expectations concerning the 
behavior of mesons in medium have been based on calculations carried out in this 
model [18-26,32]. 



2 Current-algebraic sum rules 

In this section we present a set of sum rules that are going to be explored in this paper. 
The method follows Refs. [28,30]. The sum rules follow from the SU(2)®SU(2) chiral 
charge algebra [33,34] of QCD and involve no extra assumptions, therefore are very 
general. In the context of effective chiral models such relations were derived in Refs. 
[18]. In this section we also derive the corresponding relations involving the vector 
current, i.e. involving the excitations with quantum numbers of the 5 (or ao(975)) 
meson, with I G (J PC ) = 1~(0 ++ ). For the simplicity of notation the derivation is 
made for two flavors, generically denoted by u and d. The cases involving strangeness 
(K and Kq excitations) can be obtained from the results below by replacing u or d 
by s. 

2. 1 Operator identities 

Consider the charges corresponding to vector and axial vector rotations, defined in 
the usual way as Q a = J d 3 xJ$ and Q% = $ d 3 xJ§ , with the appropriate currents 
defined as = ip^n^tp and Jg = ^7^75 ■yV'- The charges satisfy the SU(2)®SU(2) 
chiral charge algebra: 



The density of the QCD Hamiltonian is denoted as 7Yqcd- We need explicitly the 
mass term, 7i mass = ipAAtp, where the current mass matrix is M. = diag(m u , rrid). 
The canonical anticommutation rules for the quark operators, {ip a (x,t),ipJy,t)} = 
5 3 (x — y)5 al 3, and the explicit form of H mass result in the following operator identities: 



[Q», Q b ] = ie abc Q c , 



[QlQ\]=ie abc Q c . 



(1) 



[Q\ [Q 6 ,^qcd(0)]] = #))[t72, [r b /2,M]}m, 
[Ql [Ql Wqcd(0)]] = ^(0){r°/2, {r b /2,M}}m- 



(2) 
(3) 



Rewriting these relations for the neutral and charged channels (with r 
we obtain the following operator identities: 



± 



(r^r 2 )/^), 



3 



[Ql [Ql Wqcd(O)]] = m u uu(0) + m d dd(0), (4) 

[Qs, iQt,n QC B(0)]] = \{m u + m d ) (tz«(0) + dd(0)) , (5) 

[Q°,[Q°,^qcd(0)]]=0, (6) 

[Q~, IQ + , Wqcd(O)]] = \{m u - m d ) (uu(Q) - dd(0)) . (7) 

Relation (6) is trivial, since the third component of isospin is a good symmetry even 
when m u ^ m d . Relation (7) is nontrivial only if m u ^ m d . 



2.2 Gell-Mann-Oakes-Renner relations in medium 



In Appendix A we present a detailed derivation of sum rules from the above operator 
identities, and the reader is referred there for the details. The sum rules are obtained 
by the usual technique: Identities (4-5, 7) are sandwiched by a state |C), given below. 
Then, a complete set of intermediate states is inserted in the LHS of the identities. 

The state \C) is chosen to be a uniform, translationally-invariant state describing 
nuclear matter. It has fixed baryon number density, ps, and isospin density, pi=\. 
We choose to work in the rest frame of nuclear matter. Let us explain the notation 
used below: states \j a ), where a labels isospin, denote all states that can be reached 
from the state |C) by the action of the appropriate current. For instance, in the 
case of the J® operator, the states have quantum numbers of the neutral pion, 
and include all possible modes excited "on top" of nuclear matter: the vacuum pion 
branch, collective modes, lp-lh, 2p-2h, etc., excitations of the Fermi sea, etc. As 
shown in Appendix A, the sum rules involve intermediate states with momentum 
in the nuclear matter rest frame. The quantity Eja denotes the excitation energy of 
the state \j a ) (in the rest frame of nuclear matter). The symbol includes the sum 
over discrete states, as well as the integration over continuum states. 

Relations (4-5) result in sum rules which are in-medium generalizations of the Gell- 
Mann-Oakes-Renner (GMOR) relations [35]: 



-m u {uu) c -m d (dd)c = Y,Yf\ J lM\C) 

3° 

-{m u + m d )(uu + dd) c = J2 sgn(£7,-) | <J ~ I (°) I C> 

+ J2sgn(E j+ ) (j + |J+ (0)|C) 



(8) 



(9) 



Indeed, in the case of the vacuum, |C) = |vac), we can single out the one-pion 
contribution in the RHS of Eq. (8-9). Let us denote this state (with 3- momentum 0) 
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as \n a ). For example, in the case of (8) we then find (tt°\ J® (0)|vac) = m^F^o, where 
m n o and F n o are the neutral pion mass and decay constant. Therefore, we can write 

- m u (uu) vac - m d (dd) vac = mloFfr + £ |(/| J 5 ° (0)|C)| 2 = rrfaF^ + O(m 2 ).(10) 



The symbol X^o^ denotes the sum over all contributions other than the one-pion- 
state, e.g. three pions, p+n, etc. As is well known, such contributions are chirally sup- 
pressed [36]. They are also infinite, hence require renormalization. Note, however, that 
since no extra divergencies are introduced by nuclear matter, the vacuum-subtracted 
sum rules (8-9) are well-defined: 



-m u ((uu)c - (uu) vac ) - m d ([dd) c - (dd) va ^ 
= E |0'Ko(0)|C> 2 - E |0'Ko(0)|vac> 



(11) 



j" 3" 

and similarly for the sum rule (9). 
2.3 Additional sum rules 

Repeating the steps of the previous section on Eq. (7) we arrive at the sum rule 



— (m u — rrid)(uu — dd)c — vac 
= E88 n (^-) <rl^(0)|C> 2 + Esgn(^ + ) 0' + l^o + (0)|C) 



— vac, 



(12) 



where vac means the vacuum subtraction as in Eq. (11). Here the intermediate states 
have quantum numbers of the S meson, (J G (J PC ) = 1 _ (0 ++ )). 

Subsequent sum rules are obtained from Eqs. (1). The derivation repeats the steps of 
App. A. We obtain two sum rules involving the isovector density, pi=\ = ^{u^u—d^d)c'- 



(^u-Sd) c =Y,T^ <m7o(o)|c> -E 



(„t u _ d t d)c=s (r\j -(o)\c) -E 



0' + Ko(0)|C> 
(j + \4(0)\C)^ 



(13) 
(14) 



Sum rules (13) and (14) involve excitations with the quantum numbers of ir and 
5, respectively. These sum rules require no vacuum subtraction, since the left-hand 
sides involve the matrix element of the conserved vector-isovector current ( pi = \ = 
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(C\ Jq(0)|C)). If the state \C) is isosymmetric, i.e. pi=\ = 0, then the above relations 
are trivial and just reflect the isospin symmetry of the excitation spectrum. 



3 Formal results from the sum rules 



For the discussion of this section it is convenient to rewrite the sum rules using the 
identities 



^0' a |^o(0)|C) = (j a | [tf QC D, J 5 a ,o(0)J \C) = (j a mO)^{r a /2,M}m\C), 
E ja (f\JS(0)\C) = (j a \ [H Q cB,JS(0)} \C) = (f\m [r a /2,M]m\C). (15) 

Then Eq. (8) becomes 

1 2 

- m u (uu) c - m d (dd) c - vac = ^ [m u ui^u - m d di>y 5 d) \C) - vac, (16) 

3° ^l 



where j° labels all excitations with the quantum numbers of n°, Eqs. (9,13) give 



■(uu + dd) c — vac 



(17) 



=E 



2p/=i = E 



m u + m d 
2\E j -\E j - 

(m u + m d y 



(j \di l5 u\c) +Y. 



m u + m d 



2 \Ej+ \ Ej 



■j+\^3 + 



(j + \m l5 d\C) 



vac, 



2\E j -\ 3 

3 



(r\di^u\c) -y: 



(m u + m d ) 



2\E j+ \ 3 



(j + \ui lb d\C) , (18) 



where j* 1 label all excitations with the quantum numbers of 7r ± , and finally Eqs. 
(12,14) give 



— (uu — dd)c — vac 



=E 



m u - m d 



2\Ea- \ Ea- 

3 



m« - m d 



j+ 2 \ E j+\ E 3- 



(j + \ud\C) 



vac, 



2p/=i = E 2 |.|3 



0-|d«|C7> -E 



(m u - m d )' 



2\E, 



+ 



(j + \ud\C) 



(19) 



(20) 



3~ 3 

where ^ label all excitations with the quantum numbers of 5 ± . 

We stress that the above sum rules are valid for all values of current quark masses, 
i.e. not necessarily in the chiral {m u + m d — > 0) or isovector {m u — m d — * 0) limits, 
and hold for all densities ps and pi=\. 
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3. 1 Chiral limit at finite density in is o spin- symmetric medium 



Now we are going to explore several formal predictions following from Eqs. (16,17,19). 
The method has been discussed in detail in Refs. [28,30]. First, we analyze the case 
when the state |C) carries no isovector density, such as the vacuum or symmetric 
nuclear matter. To simplify notation we take the strict isovector limit m u = = m. 
In this case of exact isospin symmetry of the Hamiltionian, as well as of the state 
|C), the excitation spectrum is invariant under isospin rotations, and clearly Ejo = 
Ej+ = Ej-. Also, (uu)c = (dd) c = (qq)c- Sum rule (16) becomes 

- 2(gg)c - vac = H 0'°l (ui^u - di^d) \C) 2 - vac. (21) 

j0 E j0 



As long as the chiral symmetry is broken, (qq)c is non-zero in the chiral limit m — > 0. 
As already mentioned, the vacuum subtraction terms are of order m, thus are chirally 
small. Therefore, to match the chiral dimensions on both sides of Eq. (21), there must 



exist a state, denoted as 7r°, for which 



m 



[ir°\ (ui^u — di^^dj \C) 



1. Since the 



matrix element (ir°\ (ui'j^u — di^d) \C) is finite in the chiral limit, then it follows 



that E n o ~ vrn. Thus, we obtain the same chiral scaling as in the vacuum, where in 
the chiral limit we have = ^2fn(qq)o/ F n ~ \fw. By isospin symmetry we have 



E w o 



Ej, 



'm. 



(22) 



Note that this result is true for any baryon density. 

In principle, in the dense medium there could be more than one state contributing to 
the sum rule (21) in the chiral limit. It is known that many-body effects of the Fermi 
sea can induce additional branches of excitations, and we could have several states 
scaling as (22). Whether or not this occurs is a complicated dynamical issue. The 
formal result states that there exists at least one state scaling as (22) in the chiral 
limit. 



3.2 Chiral limit at finite density in isospin- asymmetric medium 



As shown in Ref. [28], in medium with finite isovector density, p/ = i ^ 0, the behavior 
of charged excitations in the chiral limit is radically different from the isosymmetric 
case (22). First, an obvious remark is that since the medium state |C) breaks the 
isospin invariance, the isospin symmetry of excitations is broken. In fact, at low 
densities [18] one can relate the splitting of E v+ and E n - to the Weinberg- Tomosawa 
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term in the tt-N scattering, and obtain 



E n+ - E n _ - (23) 

"it 



In this approach one takes the low-density limit prior to the chiral limit. Eq. (23) 
shows that for negative p I=1 at small densities we have E w _ > E w+ . However, Eq. 
(23) cannot be used at large densities. 

Now, following Ref. [28] , we assume that the isospin density is fixed, and employ sum 
rules (17-18). Without loss of generality we can assume that pj = \ < 0, as in the case 
of neutron stars or large nuclei. Since pi=\ is an external property of the system, i.e. 
independent of the chiral parameter, it is treated as large (finite) in the chiral limit. 
Then, also the isovector chemical potential, pi=i, defined as the minimum energy 
needed to lower the isospin by one unit, is finite in the chiral limit. 1 The excitation 
energies of positive (negative) isospin can now be written as Ej± = pi=\ ± SEj±, 
with 5Ej± > by definition of the chemical potential. Hence, for the medium with 
Pi=i < 0, or p,i = \ < 0, only the positive-isospin excitation energies can vanish. The 
negative-isospin excitation energies cannot vanish, including the case of the chiral 
limit. Hence we arrive at the result that for pi=\ < the negative- charge excitations 
are chirally large, 

Ej_ ~ 1. (24) 



Since the matrix elements (j a \ui'y 5 d\C) are non-singular in the chiral limit, the sum 
over negative-isospin excitations, j~, in Eqs. (17-18) does not contribute as m u + 
nid ~^ 0. Therefore in the chiral limit the negative-charge states do not contribute 
at all to the sum rules (17-18). For the sum rules to hold, there must exist positive- 
charge states scaling appropriately in the chiral limit [28] . Assuming there is only one 
such state, labeled n + , we find that in the chiral limit 



- (uu + dd) c = ™; + ™ d (n + \m l5 d\C) * (25) 



1 Although we have not been able to prove this statement from first principles, one can 
present a number of physical arguments in its favor. In the /^-exchange model discussed 
in Ref. [28], when an object of isospin ^3 is placed in the isospin-asymmetric medium, 
the energy gain is equal to g 2 p /m 2 p pi = \I^, and the corresponding chemical potential is 
M=\ = g 2 p /m 2 p pi=\ ~ pi=i/(2F%) , where the last equality follows from the KSFR relation. 
This shows that finite pi=\ in the chiral limit implies finite m=\. Another example is 
provided by the Fermi-gas model discussed in this paper. The expression for the chemical 
potential is (see Sec. 4 for notation) pi=\ = p + + M'£— \Jk"\ + Mj , and it is finite 
when k u / kj, i.e. when pi = \ is finite. 
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-2p/=i = 



(m u + m d ) 
2|£ 7r+ | 3 



{-n + \ui^zd\C) 



(26) 



Equations (25-26) give immediately 



{n + \ui^^d\C) 



(m u + m d ) (uu + dd)c 
2pJZ 

(m u + m d ) (uu + dd) c 



3 



0{m u + m d ), 



0(m u + m d ). 



(27) 



(28) 



This is totally different from the "usual" behavior in the chiral limit, Eq. (22): the 
excitation energy of the n + mode scales as the current quark mass itself, m u + m d , 
and not y/m u + m d . 

The formal case where more than one state contributes to the sum rules (17-18) in 
the chiral limit has been analyzed in Ref. [28] in the following way: Assume that 
the excitation energies of one of these modes scales as E n+ ~ [m u + m d ) a , and 
the corresponding matrix element scales as |(7r + |m75d|C)| ~ (m u + m^f . Since the 
matrix element is not singular in the chiral limit, one has /3 > 0. The RHS of the sum 
rule (18) contains, in the chiral limit, only the positive-charge contributions, which are 
negative-definite. Therefore no cancellations between contributions of various modes 
can occur in the chiral limit. The requirement of matching chiral powers on both sides 
of Eq. (18) gives 3a = 2 + 2/3, from which we conclude that a > 2/3. The positive- 
charge contributions to the RHS of the sum rule (17) may contain both positive and 
negative terms, since the sign of Ej+ is not constrained. This means that there can be 
cancellations of the leading chiral powers of the positive-charge modes contributing 
to the sum rule in the chiral limit. Matching of the chiral powers on both sides of 
Eq. (17) yields therefore the inequality > 1 — 2a + 2(3, which means, together with 
the previously-derived equality 3a = 2 + 2(3, that a < 1. In case of no cancellations 
a — 1. Combining the above results one gets [28] 



Note that these inequalities are nontrivial, since in the vacuum the corresponding 
power is a = 1/2 < 2/3. 

Carrying a similar analysis for the neutral pionic excitation in isospin-asymmetric 
medium we find that in the chiral limit it scales the usual way according to Eq. (22), 
i.e. as in the isosymmetric case. This is clear, since the medium does not break the 
third component of isospin. 

We conclude this section with several comments. Firstly, the behavior of charged 
pionic modes in the chiral limit is radically different when the medium breaks the 



2/3 < a < 1. 



(29) 
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isospin symmetry. For pj = \ < the negative- charge excitation has finite energy in 
the chiral limit, Eq. (24), and the positive-charge excitation becomes very soft, with 
excitation energy scaling as Eq. (27) (if there is a single mode contributing in the 
chiral limit) or according to Eq. (29) (if there are more modes). Furthermore, as we 
will show in the model calculation in the following sections, the nature of this soft 
mode can be quite complicated: it need not be the excitation branch connected to 
the pion in the vacuum, but the spin-isospin sound mode, resulting from collective 
effects in the Fermi sea. 



3.3 Isovector limit at finite density 



Now we pass to the analysis of Eq. (19-20) in the isovector limit m u — m d — > 0. To 
our knowledge, results of this section are novel. The method is the same as in the 
previous sections, with the obvious difference that now we compare the powers of 
TTifi — ma rather than m u + m d on both sides of the sum rules (19-20). 

In isosymmetric medium (uu — dd)c ~ m u — m d . Thus, the powers of m u — m d in Eq. 
(19) are matched with 

E 5 » -O(l), a = 0,+,-. (30) 



This result is not surprising: since 5 is not a Goldstone boson, there is no reason for 
its mass to vanish in the vacuum, or in isosymmetric medium. 

In asymmetric medium the situation is different. Since pj=\ = (u^u — d^d)c is large in 
the isovector limit, i.e. does not vanish when m u — m d — > 0, also (uu — dd)c is large. 
This is a direct effect of the asymmetry of the Fermi sea, since the Fermi momenta 
of u and d quarks (or protons and neutrons) are not equal. As a result, the Fermi 
seas of particles with opposite isospin contribute differently to (uu)c and (dd)c- An 
explicit example is provided later on in Eq. (35). 

Repeating the steps of Sec. 3.2 we now find that for p I=1 < there must exist at 
least one state 5 + whose excitation energy scales as m u — m d in the isovector limit. 
If there is only one such state, then we find 



(m u -m d )(uu- dd) c , n , , 
E s + = — 0(m u -m d ), 



2p 



i=i 



(5 + \ud\C) 



(m u - m d ) (uu - dd) c 



2p?=i 



0(m u - m d ), 



(31) 
(32) 



Equation (31) shows that there is an exact zero mode in the case of the strict isovector 
limit m u = m d . Such collective modes have been known to occur in many-body physics 
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[37-39]. If more such states exist, then, in analogy to the case of pionic excitations 
(cf. the discussion above Eq. (29)), we find that the excitation energies of these states 
scale as (m u — rrid) a , with 2/3 < a < 1. This result is nontrivial, since in the vacuum 
the corresponding power iso; = 0<2/3. 



3.4 Comments 

We end the formal part of this paper with several comments. We stress again that 
the sum rules of Sec. 3 are valid for arbitrary values of the current quark masses, 
not necessarily in the chiral or isovector limit, and for arbitrary densities. All kinds 
of intermediate states contribute to the sum rules: quasiparticles (poles), which can 
come in multiple branches, 2p-2h continuum, etc. 

Another remark concerns the sign of the excitation energy of a mode. As already 
noticed in Refs. [28,30], the charged excitation may have negative excitation energy. 
Note, however, that this does not mean that the system is unstable. This is because 
charged excitations change the isospin of the system. Suppose we request the state 
|C) to be the ground state of matter with isospin constrained to the value I3. A 
charged excitation in the sum rules involves isospin Is±l. Thus, its isospin is outside 
the constrained value, and even if the energy of the mode is lower than the energy 
of |C), it does not mean instability of |C). As a corollary, we notice that states of 
negative excitation energy cannot be the soft modes of Eq. (27) in the chiral limit. 
We conclude this from Eq. (27). Since the quark condensate is negative, we get (for 
media with negative isospin density) that E n + > 0. 



4 The Nambu-Jona-Lasinio model 

In the remaining parts of this paper we are going to illustrate in detail the general 
results discussed above by a model calculation. We will consider quark matter in 
the Nambu-Jona-Lasinio model [31]. This model acquired great popularity in recent 
years as a framework for calculations of meson and baryon properties, also in the 
nuclear medium [18-26] described as a Fermi gas of quarks. Although the description 
of the nuclear medium by a Fermi gas of quarks is certainly not realistic (unless, 
perhaps, at very high densities where nucleons deconfine), the model is well-suited 
for our theoretical purpose: it is consistent with the constraints of chiral symmetry. 
Therefore it complies to the current-algebra relations, and the sum rules of Sec. 3 
hold. 

The Lagrangian of the SU(2) ® SU(2) Nambu-Jona-Lasinio model is 
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C = qii^d, -M)q+ ( f {{qqf + (qh 5 r a q) 2 ) + ^ ((qr a q) 2 + (gi 75 g) 2 ) 

{(Ql^qf + (g 7 57^g) 2 ) - ^f{q lp qf , (33) 

where q is the quark field, Ai is the quark mass matrix, and G"s denote the coupling 
constants in various channels. Using the usual technique of the Hartree approximation 
one arrives at self-consistency equations for the values of the scalar-isoscalar field S, 
the neutral component of the scalar-isovector field 5, the time component of the 
neutral vector-isovector field p, and the time component of the vector-isoscalar field 
uj: 



m u + m d - m u -m d 
b = G a {uu + dd), d = 

p = 2G P (u + u — d + d), uj — Guj {u + u + d + d) . 



— Gs (uu — dd) , 



(34) 



For the numerical study in the examples below we use the following two parameter 
sets, fixed by meson properties in the vacuum: 



G (T = 7.55GeV~ 2 , G s = 5.41GeV- 2 , G p = 7.09GeV 



-2 



A = 750MeV, m u = 2.52MeV, m d = 4.52MeV, (I) 

G CT = 4.35GeV- 2 , G 5 = 3.34GeV~ 2 , G p = 12.4GeV" 2 , 
A = 954MeV, m u = 1.03MeV, m d = 3.03MeV. (II) 

The value of G^ and the uj field are not relevant, since we will look for excitations 
carrying no baryon number. Parameter set I has been used in Ref. [40] to fit the 
mesonic properties: m n , F n , m p = 765MeVQ and m v = 519MeV. This fits four 
parameters out of original six. The remaining two parameters are chosen in such a 
way that S = 361MeV, and the current quark masses m d and m u are arbitrarily split 
by 2MeV. Parameter set II also fits m^, F w , m v and S = 361MeV, but not ra p . It 
uses a much larger value for G p . Such larger values are needed if one wishes to fit the 
an scattering length [41]. Following Ref. [40], we regularize the model using the 
sharp 3-momentum cut-off. Our results do not qualitatively depend on the choice of 
the regulator, since the constraints of current algebra are satisfied. The 3-momentum 
cut-off obeys these requirements, in particular it leads to correct Ward identities [18]. 

The scalar and vector densities of the u and d quarks are equal to 



2 Note that in this fit the p meson lies just above the qq production threshold, hence the fit 
is somewhat problematic. However, this issue is not of much relevance for our illustrative 
application of the model. 
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/rl° k M 
, 9 r 3 j= = = (®(k u ~ \k\)- 0(A - |*|)) , 

(dd) = 2N C / d (e(fc d - - 0(A - |fc|)) , 

(u+u) = 2N c J-0^Q(k u -\k\), (d+d) = 2N c J^ 3 Q(k d -\k\), (35) 

where A is the sharp three-momentum cut-off, and k u and kd are the u and d quark 
Fermi momenta, and is the step function. We have introduced scalar self-energies 
of u and d quarks, given by 

M u = S + 5, M d = S-5. (36) 

Self-consistency requires that the quark propagators be evaluated with mean-fields 
(34): 

S u/d =P~ ^o(±t| + - M u/d + is sgn(fi u/d - p ) , (37) 
where jj, u and jj, d are the chemical potentials of the u and d quarks. 



5 Mean fields in medium 



We introduce the x and y variables, 



Pb Pu + Pd Pd / QO x 

x — — = , y = , (38) 

Po N c p p u + Pd 



where po — 0.17fm~ 3 is the nuclear saturation density, ps is the baryon number 
density, and p u / d are the quark densities. The variable y measures the isospin asym- 
metry of the medium. In symmetric medium y — |, and in pure neutron mat- 
ter p d = 2p u and y — |. The isospin density of the system can be written as 
pi=i = \{p u — Pd) = N c p x(j — y)- In our study we fix the x and y variables, hence we 
examine properties of quark matter at a given baryon density and isospin asymmetry. 

The first task is to find the mean fields by solving Eqs. (34). For the field p, which 
couples to the isospin current, we get immediately p = AG p p\ = i = 4G p iV c pox(| — y). 
The values of S and 5 are determined by solving numerically the first two of equations 
(34). Results for M u , M d and —\p are displayed in Fig. 1. For isosymmetric matter 
(top row) M u is practically equal to M d , and the small splitting is caused by the 
current quark mass difference, m d — m u = 2MeV. For y = | (middle row) and x in 
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Fig. 1. The quark scalar self-energies, (solid line) and M u (dashed line), and the vector 
self-energy of the d quark, —\p (dotted line) plotted as functions of x = Pb/po for several 
values of y. The top plots are for y = 1/2 (isosymmetric matter), the middle plots are for 
y = 2/3 (pure neutron matter), and the bottom plots are for y = 1 (pure d-quark matter). 
The plots on the left (right) are for parameter set I (II). 

the range 1 to 4 we find that M u is greater than by 10 — 20MeV. At maximum 
asymmetry (bottom row) the u quark is heavier than the d quark by 100 — 200MeV. 
There is a simple physical argument why M u > at y above |: the d quarks 
are more abundant, and it is energetically preferable for the system to make them 
lighter .P] 

We note that the field S has a large value of the order 0(1) if the chiral symmetry 
is broken. Otherwise it is of the order 0(m u + m^). The field 5 is large only if 
(uu — dd) ~ 0(1), which occurs in isospin-asymmetric medium. In isosymmetric 

3 Note that although the values of constituent quark masses in the two lower rows of Fig. 1 
decrease with density, in the strict sense it does not mean chiral restoration. This is because 
chiral symmetry cannot be restored when isospin is broken. Indeed, if Qz\C) ^ 0, then by 
charge algebra [Q$, Q^]\C) ^ 0, hence we cannot restore chiral symmetry, in which case 
we would have Q%\C) = 0, a=l,2,3. 
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medium (uu — dd) ~ 0{m u — m^), and 5 is small, of the order 0{m u — rrid). 



6 Meson propagators in medium 

As explained in Appendix A, only excitations "at rest" enter the sum rules. Further- 
more, we shall only consider the charged meson propagators, since the interesting 
effects take place for that case. 

In the case of no vector-isovector interactions {i.e. G p — 0), the one-quark-loop inverse 
pion propagator acquires a simple form 1 — G C J^^ where 

r d 4 k 1 1 

In presence of vector-isovector interactions there is a complication due to the well- 
known mechanism of mixing of 7r and the longitudinal component of the A\ meson. In 
that case in order to find excitation energies one has to find zeros of the determinant 
of the inverse n — A\ propagator matrix, D w (see e.g. Ref. [18] for details concerning 
this problem). The explicit form of the determinant is given in Eq.(B.4). 

It is worthwhile to look at the analytic structure of D n , or equivalently, J nn , in the 
variable q . The matter state |C) consists of the Fermi seas of d and u quarks, with 
kd > k u , as well as of the Dirac sea occupied down to the cut-off A. A positive-charge 
Fermi sea excitation moves a quark from the occupied d level to an unoccupied u 
level. Pauli blocking allows this when 

p+ ^/kj + Ml- y/k* + Ml < q < p+ y/k* + Ml- ^kl + Mj. (40) 

Thus, within these boundaries D 7T (q ) possesses a cut. The cuts associated with the 
Dirac sea are within the boundaries: 

p -jA2 + M*-y/A* + lW<to<p-y/kl + M*-y/kl + MZ, 
p + y /kl + M* + y/kl + M2<q <p+ VA 2 + M2 + + MJ. (41) 

In the delta channel we proceed analogously. We define 

Jss(q) = -iTr J ^yM k + ^W* - \<l)- (42) 

For the case G p = the inverse charged 5-meson propagator is 1 — GsJss- For finite 
G p there occurs mixing between the S meson and the longitudinal component of the 
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Fig. 2. Charged pion in isosymmetric matter (y = 1/2), parameter set I. (a) Charged pion 
excitation energy plotted as fuction of x. (b) The contribution of the pion poles to sum rule 
(9): the solid line shows the contribution of tt + (which equals to the contribution of ir~), 
the dash-dotted line shows the combined contribution of the two poles. 

p meson. This mixing is proportional to the mean field S, hence it is small, of the 
order of 0(m u — m^) in isosymmetric medium. The stated behavior can be promptly 
seen from Eq. (B.8). If the medium is asymmetric, then the mean field 5 is large, and 
such is the 5 — p mixing. The explicit form of the appropriate determinant, D$, is 
given in Eq.(B.9). The location of the cuts of D$ is of course the same as in the pion 
case. 



7 Mesons in symmetric matter 



Figure 1 shows the results of the numerical calculation of the charged pion excitation 
in symmetric matter. In Fig. 2(a) we show the position of the charged pion excitation 
at rest (usually called the in-medium pion mass) as a function of baryon density. 
The behavior is the expected one [22] , with the pion mass increasing slowly with the 
baryon density up to about x ~ 2. Above this point chiral symmetry is restored, i.e. 
S ~ 0{m u + rrid), (cf. the upper left Fig. 1) and the pion mass grows more rapidly. 

Figure 2(b) shows the anatomy of the in-medium GMOR sum rule (9). We note that 
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Fig. 3. Charged 5 meson in isosymmetric matter (y = 1/2). (a) The excitation energy 
plotted as fuction of x. (b) The contribution of the 8 poles to sum rule (12): the solid line 
shows the contribution of 5 + (which equals to the contribution of <5~, the dash-dotted line 
shows the combined contribution of the two poles. 

up to x = 2 practically all of the sum rule is saturated by the charged pion poles. 
At larger x some small ( a few per cent) strength is carried by the cuts (cf. Eqs. 
(40-41)). We have verified for all other cases shown in this paper that the sum of all 
pole and cut contributions to the sum rules adds up to 100%. This serves as a check 
of the numerical calculations. 

The case of the charged 5 excitation is displayed in Fig. 3. This excitation emerges as 
a bound state from the qq continuum at x ~ 0.6. Its mass decreases with the baryon 
density up to x ~ 2.5, and then starts growing (Fig. 3(a)). The contributions to the 
sum rule (19) are shown in Fig. 3(b). We can see, especially at lower values of x, that 
the pole contribution fall short of saturating the sum rule. Continuum contributions 
carry about 50% at x — 1 and about 15% at x > 3. 



8 Mesons in asymmetric matter 



In this section we come to the central part of our paper. We will show that in our 
model the sum rules from Sec. 3 are, for the case of isospin asymmetric medium, 
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Fig. 4. The real part of the pion determinant for Gs = 4.7GeV 2 (dashed line) and 
G 5 = 9.4GeV~ 2 (solid line), plotted function of q° for y = |. Other parameters 

are: G a = 11.7GeV~ 2 , A = 619MeV, and m u = m d = 5.7MeV. For large values of G s a 
zero of the inverse pion proparator is induced in the vicinity of the du, denoted by a blob. 
This indicates the presence of the "spin-isospin sound" . 

satisfied in a non-trivial way. This involves a collective state, specific for asymmetric 
medium. As explained e.g. in [42,43] in the framework of conventional nuclear physics, 
it is possible for the pion propagator in neutron matter to have an additional pole at 
very low excitation energies. Such an excitation is known as the spin-isospin sound. 
We will show that this phenomenon occurs in our model. 

The existence of collective modes in our model is related to the presence of the Fermi 
sea cut (40). Figure 4 shows the real part of the determinant D w for y = | (pure 
neutron matter) and two sample parameter choices, plotted as a function of the 
energy variable q° in the region of the Fermi sea cut. Let us first look at the solid 
line, corresponding to parameters with a large coupling constant Gs- The presence 
of the cut manifests itself by the two cusps. The imaginary part of D n is nonzero 
in the region between the two cusps, and vanishes outside. We notice that a zero of 
D n exists in the vicinity of the cut, indicated in the figure by a blob. This zero, at 
q° = 23MeV, corresponds to the energy of the spin-isospin sound mode. The dashed 
line, corresponding to lower Gs, also has cusps, but no zero of D n exists. This can be 
understood as follows: the cut region is wider and the function at the cusps acquires 
higher and lower values as the splitting of the scalar self-energies M u and Md is larger 
(c/. Eq.(40)). This splitting is proportional to the mean field 5, which increases with 
Gs, and with asymmetry y. Thus we have a critical behavior: above some critical 
values of asymmetry y and coupling Gs the spin-isospin mode emerges. We denote 
it by 7r s . In the example shown in Fig. 4 the excitation energy of n s is lower than 
the left boundary of the cut. We find this is the case for small values of the vector- 
isovector coupling constant G p . At sufficiently large values of G p the collective mode 
emerges at energies larger than the right boundary of the Fermi cut. In any case, 
the collective state lies very close to the Fermi sea cut, with excitation energy of the 
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order of lOMeV. 



In addition to the collective mode tts, there exist the usual two charged pion branches, 
7r + and 7r~, with excitation energies of the order of m % . These branches connect to the 
vacuum pion as the baryon density is lowered. Thus, depending on parameters and 
the value of y, we have, in our model, 2 or 3 branches of the charged pion excitations. 

For the charged 5 channel the situation is similar: for appropriate parameters and 
y > \, a collective mode 5s appears in addition to the usual 5 + and S~ modes. 

9 Sum rules in asymmetric medium 

In this section we show the results of our numerical study. For the case of pionic 
excitations these results have been already reported in Ref. [30] (for the slightly 
different parameter cases with m u = ). Figure 5 shows the results for the 7r 
channel for the parameter set I. Figure 5(a) shows the excitation energies of the 
usual branches, ir + and 7r~, and Fig. 5 (b) shows the excitation energy of the collective 
7i s mode. The dashed-dotted lines show the boundaries of the Fermi-sea cut, (40). 
The collective mode emerges from the cut at a low value of the baryon density. Its 
excitation energy is positive for x between 0.6 and 3.4, and negative otherwise. In 
Fig. 5(c) we show the relative contributions from the poles to the in- medium GMOR 
sum rule, Eq. (9), and the total contribution from the three poles, indicated by the 
dash-dotted line. The poles practically saturate the sum rule, leaving 1-2% for the 
cuts at large values of x. The contribution of tt,s to the sum rule (9) is of the order 
of a few per cent. Its sign follows the sign of the excitation energy in Fig. 5(c), as 
is apparent from Eq. (9). Figure 5(d) shows the relative contribution of the poles 
to the sum rule (13), and the total pole contribution, indicated by the dash-dotted 
line. We note that this sum rule is saturated by the pole at the 99.9% level - the 
cut contributions turn out to be very small. At larger value of x the collective 7Ts 
mode dominates over the other modes, and for x > 3 it practically saturates the sum 
rule. We note that the sign of the contributions is associated with the charge of the 
excitation, as is clear from Eq. (13). 

Figure 6 shows the result of a formal study of the chiral limit, m u + — > 0. For 
fixed values of y — | and x = 2 we lower the value of 



where here the superscript phys denotes the values from the parameter set (I). Figure 
6(a) shows that as the value of a is decreased, the excitation of the modes go down. 
The excitation energies of ir + and n~ modes go to finite values at a — > 0, and the 
excitation energy of n s goes to 0. Hence ir s is the chiral soft mode of Eq. (27). Figures 




m u + m d 



(43) 



a = 
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Fig. 5. Properties of charged pion excitations for y = | and parameter set I, plotted as 
a function of baryon density, (a) Excitation energies of tt + and 7r~. (b) Excitation energy 
of tts (dotted line) and the boundaries of the Fermi sea cut (dot-dashed line), (c) The 
relative contribution of 7r + , 7r~ and tts to the in- medium GMOR sum rule (9), and the 
total contribution from the three poles (dot-dashed line), (d) The relative contribution of 
7r + , 7T~ and tts to the sum rule (13). The total contribution from the three poles (dot-dashed 
line) practically saturates the sum rule. 
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Fig. 6. Study of the chiral limit for x = 2 and y = |. Convention for line the same as in 
Fig. 5. (a) Excitation energies, (b) relative contributions to the sum rule (9), and (c) relative 
contributions to the sum rule (13), plotted functions of a, Eq. (43). The spin-isospin sound 
mode its is the chiral soft mode of Eq. (27). Please, note that log a means In a. 

6(b-c)show that in the chiral limit of a — > the us saturates the sum rules (9,13). 
However, this happens at very low values of a, around In a = —4 or —5. Such values 
of a would correspond to the vacuum value of the pion mass of the order of 15MeV. 
This indicates that from the point of view of the sum rules we are quite far away from 
the chiral limit with the physical values of current quark masses, i.e. with a — 1. 

In order to better illustrate this point we show in Fig. 7 how the excitation ener- 
gies of various modes approach the chiral limit. Following Ref. [28] let us introduce 
dim(X) = linim^oQog^/ logm), which we call the chiral dimension of quantity X. 
In the chiral limit a quantity X has some scaling with a power of m u + m^. The 
function dim(X) extracts this power (for instance in the vacuum dim(m 7r ) = 1/2). 
The dotted line in Fig. 7 shows the chiral dimension of the excitation energy of tt s , 
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Fig. 7. Chiral dimensions of pionic modes plotted as functions of log a, cf. Eq. (43). For 
the symmetric matter (middle solid line, labeled y = ^, the chiral dimension of the pionic 
excitations tends to \ in the chiral limit (log a — > — oo). For asymmetric matter (y = |) 
the chiral dimension of 7rg tends to 1, and the chiral dimensions of tt + and ir~ tend to in 
the chiral limit. Please, note that log a means In a. 

which tends to 1 in the chiral limit, according to Eq. (27). The chiral dimensions of 
7r + and 7r~ go to in the chiral limit. The solid line in the middle of the plot is for 
7r + or 7r~ in symmetric matter, y — |. In that case, according to Eq. (22), the chiral 
dimension goes to | in the chiral limit. 

Figures 8-9 show the results for the parameter set II, also for y = | as in the case 
discussed above. The qualitative difference between the present and the former case 
is that now the spin-isospin sound mode lis has negative excitation energy for all 
values of x. Therefore in the chiral limit it is the ir + mode, not ns, which becomes 
the chiral soft mode of Eq. (27) (see Fig. 9 ). 

Now we pass to the discussion of the 5 channel, which is done for the parameter set 
I only, and for y — | . Figure 10 (a) shows the excitation energies of the 5 + and 5~ 
branches. After emerging from the qq continuum their energies first decrease until 
x ~ 3, and then start increasing. The collective mode 5s emerges from the cut at 
x ~ 0.75 (Fig. 11(b)). Its excitation energy is negative and small, less than 3MeV. 
Figures 10(c-d) show the relative contributions to the sum rules (12) and (14). We 
note that the 5s mode plays a major role in sum rule (12), and completely dominates 
sum rule (14). 

Figure 11 shows the isovector limit for the 5 channel. In this case 

= m u - m d 

phys phys ' V / 

ml ■' — m d 



We can see that the 5s mode is the isovector soft mode of Eq. (31). Its excitation 
energy drops linearly to as a is decreased (Fig. 11(a)), and the sum rules are 
completely saturated by the 5s mode in the isovector limit of a — > 0. 
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Fig. 8. Same as Fig. 5 for parameter set II 
10 Concluding remarks 



There are several messages which follow from our calculation. Firstly, we note that in 
order to satisfy the current-algebraic sum rules it is necessary to include all modes, in 
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Fig. 9. Same as Fig. 6 for parameter set II. The 7r + mode is the chiral soft mode of Eq. 
(27). 

particular the spin-isospin sounds. Certainly, a nuclear system is a very complicated 
object, and even our simple model, treated at the lp-lh level, has revealed a rich 
structure of the excitation spectrum. The power of the current-algebraic sum rules 
relies in the fact that they relate in a non-trivial way the properties of these excitations 
to the quark condensate and the isospin density. 

One may ask the following general questions: How far are we in a nuclear system from 
the strict chiral limit ( m u + 0) and the strict isovector limit ( m u — ma — > 0) 

in the real world, i.e. in a dense nuclear system, and with the physical values of 
m u and m^. The results shown in Figs. 7 and 11 indicate, that in moderately-dense 
isospin-asymmetric systems we are far away from the chiral limit, and very close to 
the isovector limit. From Fig. 7 we find that the its mode excitation energy scales 
linearly with m u + starting from m n ~ V e~ 6 139MeV ~ 7MeV, much lower than 
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Fig. 10. Properties of charged 5 excitations for y = | and parameter set I, plotted as a 
function of baryon density, (a) Excitation energies of 5 + and 5~ . (b) Excitation energy of 
5s (dotted line) and the boundaries of the Fermi sea cut (dot-dashed line), (c) The relative 
contribution of 5 + , 5~ and 5s to the sum rule (12), and the total contribution from the three 
poles (dot-dashed line), (d) The relative contribution of 5 + , 5~ (indistinguishable from 0) 
and 5s to the sum rule (14). The contribution of 5s saturates the sum rule. 
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Fig. 11. Study of the strict isovector limit for x = 2 and y = |. Convention for lines is the 
same as in Fig. 10. (a) Excitation energies, (b) relative contributions to the sum rule (12), 
and (c) relative contributions to the sum rule (14), plotted as functions of a, Eq. (44). The 
spin-isospin sound mode 5$ is the isovector soft mode of Eq. (31). It completely satisfies 
the sum rules in the strict isovector limit. Please, note that log a means In a. 

the physical value. On the other hand, Fig. 10(a) shows that the excitation energy 
of the 5s mode scales linearly with m u — already at physical value, corresponding 
to a = 0. 



Another comment is relevant for application of effective chiral Lagrangians to nuclear 
systems. In this approach one basically assumes that there is one pion quasiparticle 
in a nuclear medium, albeit with modified properties compared to the vacuum. In our 
model we find additional branches. Since they contribute largely to the sum rules, 
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they cannot be neglected. In an effective model they should be included as additional 
degrees of freedom. 

The final remark concerns strangeness. Although in this paper we have worked for 
simplicity with two flavors, extension to three flavors is straightforward. In fact, one 
can make a simple "translation" of the sum rules of Sec. 3 to the case of any flavor. 
For example, changing the d (or u) quark to s we obtain the case of charged (neutral) 
kaons. This is simply the replacement of /-spin by U or V spins. Note that nuclear 
matter is asymmetric with respect to U and V spins, therefore kaonic excitations 
on top of nuclear matter are parallel to the case of charged pionic excitation on top 
of isospin-asymmetric matter. We note that recently Refs. [44] discussed the kaonic 
excitations in the Fermi gas of quarks in the Nambu-Jona-Lasinio model. 



A Derivation of sum rules in medium 



In this appendix we explain the derivation of sum rules (8-14). Although the technique 
is very well known, we believe it is worthwhile to remind it in some greater detail 
in order to point out the differences between the derivation in the vacuum and in 
a medium. The first step in deriving the sum rules is to sandwich both sides of 
Eqs. (4-7) by the medium state |C). On the RHS this leads to a "known" quantity 
involving in-medium condensates (C\uu(0)\C) = (uu) c and (dd) C - Next, one inserts 
a complete set of intermediate states \k, j) between the current operators on the LHS. 
These states are eigenstates of the momentum operator, P, and of the Hamiltonian 
H = J d 3 xH(x). They can be labeled by additional quantum numbers, e.g. isospin. 
The medium state is also an eigenstate of P and H: 

P\C)=p c \C), H\C)=S c \k,j). (A.l) 



For matter of a large volume V the quantities pc and 8c are proportional to V . It is 
convenient to measure the momentum and the energy of intermediate states relative 
to the state \C), i.e. 

P\k,j) = (k + pc) \k,j), H\k,j) = (E 3 (k) + E c ) \k,j). (A.2) 



Quantities k and Ej(k) form a Lorentz four-vector. Thus the Lorentz-invariant mesure 
of integration is 5^-- , and the unit operator can be decomposed as follows [28]Q: 

\2Ej(k)\ 

1 = V / f k - \k,j)(k,j\. (A.3) 



4 The measure of integration is the same as for example in the case of phonon excitations 
on top of a solid. 
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We illustrate the method on Eq. (7). We rewrite the LHS, insert the unit operator 
(A. 3), express the charges by time-components of currents, shift the coordinates of 
the currents with the translation operator, and use Eq. (A. 2): 



[Q-,[Q + ,^qcd(0)]]=£/ 



d 3 k 



x 



((C\Q + \k,j)(k,j\[Q-,H(0)}\C) - (C\[Q-,H(0)]\k,j)(k,j\Q + \C)) = 

V / ^-^ — — / d 3 y [ d 3 x x 

(2n) 3 2\E,(k)\J •/ 

((C\4(y)\lj}(k,j\[J (x),H(0)]\C} - (C\[J (y^ 



d x x 



(2n) 3 2\E 3 (k) 

(cy pV % + (o) e -<* + ^ 

-(C\^ pc - 5 [Jo(0),n(0)]e- i ^°)-%j)(k^ 
d 3 k 



E 



5 3 (k) (C\J+(0)\k = 0,j)(k = 0,j|[J -(0)(0), d 3 xH(0)}\C) 



2\E 3 (k = 0)\ 

- (C\[J (0)(0),Jd 3 xH(0)}\k = 0,j)(k = 0,j\4(0)\C} 

1 



: E 



(E J (k = 0)(C\J+(0)\k = 0,j)(k = 0,j\J (0)\C} 



3 2|^-(fc = 0)| 
+ E,(k = 0)(C\J (0)\k = 0,j)(k = 0,j\4(0)\C)) = 

: E 1 0'" I Jo (0) I C) f + £ ^gn^. + | (j+ 1 J + (0) | C) 



In the last line we have decomposed the sum over indices j into the sum over positive 
and negative isospin excitations. We have introduced the short-hand notation \j + ) 
and for such excitations with relative momentum k — 0, and denoted their 
excitation energies by Ej+ and Ej-. This completes the derivation of the sum rule 
(12). With all other sum rules the steps are exactly the same as described above. 



B Charged meson propagators in medium 



The quark bubble for a meson channel is defined as 

Jnv(?) = -*Tr / -0fS u {k + \q)V'S d {k - ±q), 



(B.l) 
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where T = 70^70 [40]. We use T = 75 in the pion vertex, which allows to get rid 
of factors of % in Ward identities below. For the considered case of q = the tt-A 1 
mixing involves the time components of the axial propagator, , and the mixed 
propagator, J° Ai . The determinant of the inverse n-Ai propagator is equal to 

D^q ) = (1 - G a J^(q°))(l + G p jZ Al (q°)) + CG P (4( 9 )) 2 (B.2) 



The signs follow the convention for signs of the coupling constants in Eq. (33). The 
following Ward identities hold among the bubble functions [40]: 



(?o - p)JZaM°) = ZSJ\{q Q ) + 2(«t u - dU), 
(go - p)J\(q°) = 2SJ^(q°) + 2{uu + dd), (B.3) 

where p and S are defined in Eq.(34). These identities follow from the general re- 
quirements of chiral symmetry [40]. They can be explicitly verified to hold with our 
choice of the 3-momentum regulator. Using Eqs. (B.3, 34) we can rewrite Eq. (B.2) as 



Ar(<?o) 



m u + m d ( 2(m u + m^SGjGp _ ^ x 



Qo-P [ 
G a J nn (qo) 



2S \ g (go - p) 
2S - (m u + m d ) 
2S 



(B.4) 



This form is convenient, since it involves only one bubble function, J n7T , which has 
the explicit form 



d 3 k (P ~ go) + 25M U /Jk 2 + Ml 



I W (p - g ) 2 + 2(p - g ) V /P + M 2 + 4S5 
+ (u^d, p^-p, <$->-<$, g -> -g ) . 



(B.5) 



The zeros of D n (q ) correspond to poles of the mixed charged ir-A 1 propagator. The 
pole contributions to sum rules (9,13) are explicitly given by the expression 



sgn(£, ± ) (f I J± I C) 



m u + m d [SJ ww (q Q ) - G^^S - (m u + m d )/2)\ 



%- P 



dD 7T (q )/dq 



(B.6) 



In the S-p channel we obtain, if full analogy to Eq. (B.2-B.6), 

D S (q°) = (1 " G s Jss(q°))(l + G p J™(q )) + G s G p (Jl(q°)y. (B.7) 
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Through the use of Ward identities 



(ft - P)W) = ^JsM°) + 2(«t« - rfty, 

(ft - = 25J 5 ,(g°) + 2(uu - dd), (B.8) 

where p and 5 are defined in Eq.(34), we can rewrite Eq. (B.7) as 



D s (q ) 



ft \m u -m d | ( 2(m u -m d )5G 5 l G p A 



ft-p I 
GsJssiqo) 



25 V ft (ft ~ P) 

2£ - (m M - wij) 
25 



(B.9) 



where J$s is explicitly given by 



' d*k (p-qo) + 2SM u /Jk* + M* 



Jss{qo) = 4iV c / — —r ^= 



fcu , , v „ w , , . v „ wyv .. . M* + AS8 
+ (u^d, p^-p, 5^-5, q -> -g ) • (B.10) 

The zeros of Ds(qo) correspond to poles of the mixed charged 5-p propagator. The 
pole contributions to the sum rules (14) are explicitly obtained from the expression 



sgn(£ i± ; 



(3 ± I 4 I c) 



m d 



SJss(qo) ~ G 5 1 {5 - (m u - m d )/2) 



ft 



dD s (q )/dq 



QO=Ef 

(B.H) 
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